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We have measured resonance spectra in a superconducting microwave cavity with the shape of 
a three-dimensional generalized Bunimovich stadium billiard and analyzed their spectral fluctua- 
tion properties. The experimental length spectrum exhibits contributions from periodic orbits of 
non-generic modes and from unstable periodic orbits of the underlying classical system. It is well 
reproduced by our theoretical calculations based on the trace formula derived by Balian and Du- 
plantier for chaotic electromagnetic cavities. 
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In the last few decades, billiard systems have pro- 
vided a most appropriate model for the understanding 
and discussion of classical chaotic Hamiltonian systems 
[EJ and their quantum counterparts j^]. In two dimen- 
sions the Schrodingcr equation for quantum billiards co- 
incides with the Helmholtz equation for microwave cavi- 
ties of corresponding shape. This correspondence is the 
basis for experimental studies of quantum manifestations 
of classical chaos. Early experiments focused on spectral 
fluctuation properties in chaotic billiards HQ- About 
30 years ago Gutzwiller established a direct relationship 
between the oscillating part of the density of states and 
the properties of the classical periodic orbits by means 
of a trace formula [El. Gutzwiller's trace formula and its 
extension to systems with mixed classical dynamics || 
were verified experimentally using flat microwave cav- 
ities [Q0. Three-dimensional chaotic billiard systems 
have been scarcely studied theoretically Hll. There is 
no analogy between quantum billiards and electromag- 
netic cavities in three dimensions. Still, the latter are 
of great interest for the study of wave dynamical phe- 
nomena in chaotic systems. Spectral properties of the 
vectorial Helmholtz equation have recently been studied 
theoretically for integrable systems uO and experimen- 
tally for mixed and chaotic systems |ll -|l3| . 

We present the first experimental test of the trace for- 
mula derived by Balian and Duplantier |jl4) for chaotic 
electromagnetic resonators. The main difficulty to over- 
come consists in the construction of a microwave res- 
onator that is completely chaotic while permitting only 
a few non-generic modes. We meet these requirements 
by using a superconducting resonator with the shape of 
a desymmetrized three-dimensional stadium billiard [Ej| . 
Figure [f] shows that this billiard consists of two quarter 
cylinders with radii r% and r<x, resp. The electromagnetic 
resonator has dimensions r\ = 200.0 mm, T2 = 141.4 mm 
and is made of niobium which becomes superconduct- 
ing at temperatures below 9.2 K. This tremendously in- 



creases the resolution of the measured spectra due to a 
quality factor of up to 10 7 compared to 10 3 in normal con- 
ducting resonators. The measurements were performed 
at a temperature of 4.2 K for frequencies / up to 20 GHz. 
We show a typical spectrum in Fig. |^. The evaluation 
of four reflection and six transmission spectra yielded 



18764 resonances. According to Weyl's formula [|14 16 
for electromagnetic cavities the smooth part of the inte- 
grated resonance density, A smoot h(/), is a polynomial of 
third order in the frequency, where, in contrast to the 
corresponding quantum case, the quadratic term is ab- 
sent. Its coefficients have been obtained by a fit to N(f). 
The fluctuating part Nfi uc (f) is shown in Fig. it still 
displays smooth oscillations, which are due to the non- 
generic periodic orbits. Our cavity (see Fig. [I]) exhibits 
two types of non-generic orbits. First, there are two fam- 
ilies of three-dimensional, marginally stable "bouncing 
ball" orbits with length 2r\ and 2r2 that evolve parallel 
to the axis of the two cylinders. Second, there are tra- 
jectories inside the plane z — that are linearly stable 
with respect to deviations parallel to this plane and un- 
stable with respect to perpendicular deviations. While 
both types of orbits are of measure zero, they generate 
smooth oscillations in the staircase function and yield 
dominant peaks in the length spectrum, i.e. the absolute 
value of the Fourier transform of the fluctuating part of 
the fc-dependent density of states, where k = 27r//co is 
the wave number, cq denotes the speed of light. Such 
effects of non-generic periodic orbits have been found 
in various types of billiards j|,[l3 17 1. The experimen- 
tal length spectrum is shown in the top of Fig. |^. To 
obtain an analytical expression for the staircase function 
of the non-generic modes we combine the semiclassical 
method of Ref. 12 with the adiabatic method of Ref. 
Ipisj . The quantum adiabatic method |ll| has been ap- 
plied to the calculation of bouncing ball modes in the two- 
dimensional stadium billiard. This method corresponds 
to a Born-Oppenheimer approximation, where the fast 
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coordinate is parallel to the classical motion correspond- 
ing to the bouncing ball orbits, while the slow variable is 
transversal to it. Accordingly, we assume that the modes 
in the x— and y— directions are adiabatically decoupled 
from the modes in the z-direction and quantize the rect- 
angle with side lengths l x and l y . Then, the x— and 
y— component of the wave vector k ng of the non-generic 
modes are given as k x = nfi/l x and k y = itv/ly for in- 
tegers ji and v, where, according to the geometry of our 
billiard the z-dependence of the lengths l x and l y is 



l x (z) = ri, ly{z) = \Jr\- z 2 



ly(z) = r 2 , l x {z) = \/r\- z 1 



for < z < r-i 
for —v\ < z < 0. 



Combining this adiabatic method with the approach of 
rcf. JlJ] we express the staircase function N ng (k) for 
the non-generic modes (ng) as N ng (k) = Tr9(fc 2 — fc^ g ), 

where the trace is over the wave vector k ng of the non- 
generic modes 



N ng (k) = 



dzdk z 
2tt 
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Note that the integration over z is restricted to the inter- 
val [—7*1, r 2 ] while the fc z -integration is unrestricted. The 
integers fi and v label the modes in x— and y— directions, 
resp. Dirichlet (/i, v > 0) or von Neumann (/x, v > 0) 
boundary conditions correspond to the magnetic and 
electric non-generic modes, resp. 

For the electromagnetic cavity under consideration, 
the non-generic magnetic and electric modes decouple. 
Thus, the non-generic contribution to the staircase func- 
tion is given by the sum of the quantum mechanical ex- 
pression for Dirichlet and von Neumann boundary condi- 
tions. Figure^ shows that the smooth oscillations of the 
fluctuating part of the experimental staircase function 
are well described by our expression. Thus, the adiabatic 
method yields a very good approximation for the contri- 
butions of the non-generic modes to the staircase func- 
tion. In order to study the local fluctuation properties 
in the resonance spectrum, we rescaled the resonances to 
unit mean spacing and subtracted the non-generic contri- 
butions from the spectrum. For a time-reversal invariant, 
classically chaotic system the local fluctuation properties 
are expected to coincide with those of random matrices 
from the Gaussian orthogonal ensemble (GOE) M,E0[. 
We however find notable deviations and attribute them 
to a partial decoupling between electric and magnetic 
modes which is prominent at low frequencies. Indeed, 
our spectral statistics are in very good agreement with 
that obtained for random matrices H from an ensemble, 
that models two coupled, chaotic systems pi|- 



ft = ( He VXDVij 



Here, H e and H m are matrices from the GOE with di- 
mensions N e and N m , resp. The coupling matrix ele- 
ments Vij are Gaussian random variables with zero mean 
and unit variance, D is the mean level spacing, and A is 
the coupling strength. The random matrix model inter- 
polates between two uncoupled GOE-systems at A = 
and one GOE at A = 1. We obtain the best agreement 
between our experimental level statistics and that of ran- 
dom matrices H as defined in Eq. (|^) when we treat A as 
a fitting parameter and set N e (N m ) equal to the number 
of eigenvalues of the scalar Helmholtz equation with von 
Neumann (Dirichlet) conditions on the boundary of the 
billiard. This suggests an interpretation of model (||) in 
terms of electric and magnetic modes that are coupled 
with strength A. Figure || shows the level-spacing dis- 
tribution and the A3-statistics for a frequency range of 
5-10 GHz (A = 0.352) and 18-18.5 GHz (A = 0.658). The 
level-spacing distributions agree well with the GOE and 
with those of the model defined by Equation @. The 
A3-statistics coincides with that for the random matrix 
model @ but differs from the GOE. Note, however, that 
the GOE is approached with increasing frequency, (2- 
5 GHz: A = 0.07, 11.5-13 GHz: A = 0.546, 15-17 GHz: 
A = 0.633). We furthermore found that 2?rA = Y/D, 
where 1/D ist the level density, while T shows only a 
weak dependence on frequency for the frequency range 
from to 20 GHz. This is similar to a behavior that 
was observed in studies of isospin mixing in nuclei [ p4| . 
The partial decoupling of electric and magnetic modes at 
low frequencies was not observed in the three-dimensional 
Sinai billiard [jl3| and seems to be due to the particular 
geometry of the three-dimensional generalized stadium 
billiard. Its desymmetrized version (see Fig. Q) corre- 
sponds to two quarter cylinders, which are rotated by 7r/2 
with respect to each other. The electric and the magnetic 
modes are completely decoupled in each of these quarter 
cylinders [jlo|| . 

Let us finally turn to a semiclassical analysis of the ex- 
perimental spectrum. According to jl4j the generic con- 
tribution to the fluctuating part of the density of states 
is semiclassically given by the periodic orbit sum, i.e. the 
trace formula 



Pfluc(fc) = 



2 cos (4>p)L p /iT 
|det(l-M p )| 1 /2 



{kLp - jUp) ■ (3) 



(2) 



The factor 2 cos <pp stems from the polarization and is 
thus due to the vectorial character of the underlying wave 
equation. It vanishes for orbits with an odd number of 
reflections. The remaining quantities in Eq. (Q) are iden- 
tical to those appearing in Gutzwiller's trace formula [^). 
We numerically determined the first 381 periodic orbits, 
their lengths L p up to about 1.5 m, and their stability ma- 
trices M p by two different methods. Based on an ansatz 
for a symbolic code (see, e.g., Chap. 8 in ref. ||) we 
determined all orbits in the full (not the desymmetrized) 
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three-dimensional stadium billiard for up to eight reflec- 
tions off the boundary. For more reflections the numeri- 
cal evaluations become too time consuming. The result- 
ing list of periodic orbits was checked for completeness 
and enlarged by a numerical search of fixed points in the 
Poincare surface of section [z = Q,p z > 0). The lat- 
ter method is blind to some whispering gallery modes. 
As in the two-dimensional stadium billiard our device 
exhibits infinitely many whispering gallery modes. The 
shortest of these orbits propagate along the curved edges 
of the billiard and have lengths 2r\ + nr-i sa 0.843m, 
2r2 + ttti w 0.910m and 7r(ri + 7*2) ~ 1.07m. There are 
no significant peaks associated with these lengths (see 
Fig. [|, top). Similar cancellation effects have been ob- 
served in the two-dimensional stadium billiard |p5| . For 
the computation of the Maslov indices fi p we followed 
p6f . Figure ^ shows a comparison between the experi- 
mentally obtained length spectrum, the length spectrum 
for the sum of non-generic orbits [Eq. (13)] plus unsta- 
ble periodic orbits [Eq. and the length spectrum 
of the unstable periodic orbits alone. Diamonds mark 
those peaks of the non-generic contributions, that can- 
not be assigned to bouncing ball modes but to orbits in 
the plane z — 0. Evidently, the theoretical reconstruction 
describes the experiment rather well. One peculiarity are 
peaks that appear at about I ~ 0.77 m and I ~ 0.97 m. 
While we could not find corresponding orbits inside the 
billiard, there are two orbits in the boundary plane y = 0, 
whose lengths and stability amplitudes exactly match 
those peaks. Peaks at lengths below I — 0.28 m, which 
corresponds to the shortest possible periodic orbit, are 
the result of the (unavoidable) experimental inaccuracy. 
The good agreement between the experimental and the 
theoretical length spectrum quickly deteriorates beyond 
I 1.3 m. We recall the exponential proliferation of 
long orbits which enter the trace formula ([|). It might 
be that we have missed some periodic orbits in our nu- 
merical reconstruction of the length spectrum. Note that 
a considerable fraction of the short periodic orbits have 
1 2 cos cf>p | — 2. In a semiclassical picture the electric and 
magnetic modes decouple locally on such an orbit. This 
finding supports our interpretation of the level statistics 
in terms of a partial decoupling between electric and mag- 
netic modes. 

In summary, we have investigated wave chaotic phe- 
nomena in a superconducting three-dimensional mi- 
crowave resonator. Spectral fluctuations on short fre- 
quency scales agree well with those of random matrices 
from an ensemble modelling two coupled chaotic systems. 
We interpret this as a partial decoupling of the elec- 
tric and magnetic modes in the low-frequency domain. 
The length spectrum can be understood in terms of non- 
generic modes and in terms of unstable periodic orbits of 
the underlying classical system. 
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FIG. 1. Desymmetrized version of the three-dimensional 
generalized stadium billiard. The plane z = is shaded. 
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FIG. 2. Typical transmission spectrum for the frequency 
range 13-13.25 GHz, where the resonances are still well sepa- 
rated. 
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FIG. 3. Fluctuating part of the experimentally obtained 
staircase function (dashed line) compared to the fluctuating 
part of the staircase function for the non-generic modes (full 
line). We only show that range of the spectrum where both 
curves are distinguishable. 
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FIG. 4. Length spectra from the experimental data (top) 
compared to (middle) the sum of non-generic modes [from 
Eq. (|l])] plus the unstable periodic orbits [trace formula (^) 
of Balian and Duplantier], and the unstable periodic orbits 
(bottom). Diamonds mark those peaks that are assigned to 
periodic orbits in the plane z — 0. 




FIG. 5. Left figures: Experimental level-spacing distribu- 
tion (histogram) and A3-statistics (diamonds) for the fre- 
quency range 5-10 GHz (A = 0.352) compared to GOE (full 
line) and the predictions from the model (see Eq. (^) ) (dashed 
line). Right figures: Same as left figures but for the frequency 
range 18-18.5 GHz (A = 0.658). 
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